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ABSTRACT 

We present our discovery of a narrow-line Baldwin effect, an anti-correlation 
between the equivalent width (EW) of a line and the flux of the associated con- 
tinuum, in 5-20/zm mid-infrared lines from a sample of 68 Active Galactic Nuclei 
(AGN), located at z<0.5, observed with the Infrared Spectrograph on the Spitzer 
Space Telescope. Our analysis reveals a clear anti-correlation between the EW of 
the [SIV] 10.51/im, [Nell] 12.81/im, and [Nelll] 15.56/im lines and their mid-IR 
continuum luminosities, while the Baldwin effect for [NeV] 14.32/xm is not as 
obvious. We suggest that this anti-correlation is driven by the central AGN, not 
circumnuclear star formation in the host galaxy and present a new method of an- 
alyzing this effect in mid-IR lines. We also find that the slope of the narrow-line 
Baldwin effect in the mid-infrared does not appear to steepen with increasing 
ionization potential. Examining the dependence of the EW to the Eddington 
Ratio (L/LEdd) we find no strong relationship for mid-IR lines. Our study indi- 
cates that the narrow-line mid-infrared Baldwin Effect is quite different from the 
broad-line optical/UV Baldwin effect and it is possible that the two effects are 
unrelated. The discovered anti-correlations open new possibilities in the under- 
standing the physics of the ionizing region and the continuum reprocessing by 
dust. 

Subject headings: galaxies: active — galaxies: quasars: emission lines — infrared: 
galaxies 

department of Astronomy, University of Florida, Gainesville, FL 32611 
2 Cornell University, Astronomy Department, Ithaca, NY 14853-6801 
3 University of Texas, Department of Astronomy, Austin, TX 78712 
4 University of Crete, Department of Physics, GR-71003, Heraklion, Greece. 

5 IESL/Foundation for Research and Technology - Hellas, GR-71110, Heraklion, Greece and Chercheur 
Associe, Observatoire de Paris, F-75014, Paris, France 



- 2 - 



1. Introduction 

The Baldwin Effect, first discovered by Baldwin (1977), reports the decrease of the 
equivalent width (EW) of the broad CIV1549A line with increasing UV luminosity in active 
galactic nuclei. The relationship was initially established hoping quasars could be used 
as potential standard candles in observational cosmology. Extended examination of the 
relationship over the past decades for both quasars and Seyferts demonstrated that the 
relationship is not caused by selection effects, but its cosmological use is limited due to large 
scatter (see review of Osmer & Shields 1999 and Kinney et al. 1990; Wilkes et al. 1999; 
Green, Forster & Kuraszkiewicz 2001; Croom et al. 2002; Dietrich et al. 2002; Kuraszkiewicz 
et al. 2002; Shang et al. 2003). 

Significant correlations also exist between the continuum emission and the EW in other 
UV and optical emission lines including Lya, H/3, CIV, CIII, Ly/3, OIV, OI, CII, AMI, 
CIII, Mgll, and SilV+OIV. It was also found that the slope of these relationships appear to 
increase with increasing ionization potential. In addition, an X-ray Baldwin effect has also 
been reported in Fe Ka (Iwasawa & Taniguchi 1993; Nandra et al. 1997; Page et al. 2004). 

The physical origin of this effect is still not clear. A plausible explanation is the softening 
of the ionizing continuum shape with increasing L, which would lead to weaker emission lines 
compared to the local continuum (Baskin & Laor 2004). However, this has been challenged 
by Wilkes et al. (1999), who found no correlation between any of the UV and optical lines 
with the X-ray luminosity or X-ray slope. They suggest a model in which limb darkening 
and the projected surface area of an optically thick, geometrically thin disk combine to cause 
the Baldwin effect. 

Some have also argued that the Eddington ratio L/LEdd, a tracer of AGN accretion, 
may drive the Baldwin effect (Boroson & Green 1992). The Baldwin effect may then just 
be a secondary correlation induced by the tendency of more luminous AGN to have a higher 
L/L Edd (Baskin & Laor 2004; Shang et al. 2003). Others have argued that the fundamental 
driver is the mass of the supermassive black hole M BH instead of L or L/ L Edd (Warner et 
al. 2003, 2008). One further plausible explanation is that metallicity of the gas in the AGN 
affects the EW of the lines in a way that would generate a Baldwin Effect (Deitrich et al 
2002). 

Most of the discussion of the Baldwin Effect has focused on broad-lines, but a few 
papers have also noticed a narrow-line Baldwin effect (Green et al. 2001, Croom et al. 2002, 
Boroson and Green 1992). Since the narrow-line region in quasars may extend to kpc scales, 
the physics related to the narrow-line Baldwin effect may be different from those driving the 
broad-line Baldwin effect (Osmer and Shields 1999), and may simply be due to the covering 
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factor of the narrow-line region (Page et al. 2004). This is manifest in a "disappearing 
NLR" model, where the NLR size is related to the AGN luminosity and highly luminous 
AGN would have weak or even non-existent NLR (Croom et al. 2002). 

Complications arise when moving into the mid-infrared. Several authors have noted 
that in many cases, the IR spectra of AGN do not reflect their optical or UV classifications 
as reprocessing of the ionizing radiation by the intervening dust as well as circumnuclear 
star formation activity affects the mid-IR spectral features (Lutz et al. 1998, Laurent et 
al. 2000, Armus et al. 2007, and Spoon et al. 2007). As a result, quantifying accurately 
the AGN contribution to the infrared or bolometric luminosity of dust enshrouded galaxies 
is still largely unanswered problem (see Charmandaris 2008 for a review). Ascertaining the 
extent of AGN domination in the mid-IR presents unique challenges to determining the 
extent of a Baldwin Effect. We embarked on the detailed study of the effect using Spitzer 
data (Keremedjiev & Hao 2006) while an analysis using ground based observations with the 
VLT/VISIR and comparing to x-ray luminosity has been presented by Honig et al. (2008). 

In this paper, we report on our discovery of a narrow-line Baldwin effect in the mid- 
IR based on Spitzer observations of a large sample which consists of 68 optically classified 
AGN. Our observations and data analysis are presented in Section 2, our results and detected 
correlations are shown in Section 3, while the implications of those are discussed in Section 4. 



2. Observations and Data Reduction 

We compiled a sample of 16 PG Quasars, 33 Seyfert Galaxies, and 19 2MASS AGN with 
z<0.5 observed with the Spitzer Space Telescope's Infrared Spectrograph (IRsfl All objects 
were observed in the short-high (SH), short-low (SL) and long-low (LL) IRS modules. The 
low resolution modules SL (5.2-14.5/im) and LL (14.0-38.0/im) have a spectral resolution of 
64-128, whereas the high resolution module SH (9.9-19.6/zm) has a resolution of ~ 600 (Houck 
et al. 2004). For this experiment, our effective coverage is limited by the SH bandpass. The 
data were reduced at the Spitzer Science Center (SSC) using data reduction pipeline 13.2. 
The Basic Calibrated Data (BCD) products from the SSC were used in our analysis. 

To obtain the low resolution spectra, we co-added the images of the same nod position. 
The background was then subtracted from the co-added image of the other nod position. 
The mid-IR spectra were extracted from these images with the Spectral Modeling, Analysis, 



1 The IRS was a collaborative venture between Cornell University and Ball Aerospace Corporation funded 
by NASA through the Jet Propulsion Laboratory and the Ames Research Center. 
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and Reduction Tool (SMART Ver. 5.5.1 Higdon et al. 2004) using a variable width aperture 
to recover the diffraction limited point-spread function. 

For high resolution data, the background cannot be subtracted in the same way as in 
the low resolution cases. Due to the limited size of both apertures of IRS high resolution 
modules, even point sources tend to fill them making differencing between the two nods 
impossible. Thus, images of the same nod position were co-added using median averaging and 
full aperture extractions were conducted in SMART. The resulting spectra were cosmetically 
trimmed by eye for order overlap across all orders. This was done to compensate for the 
low S/N present in most spectra. Unfortunately separate background observations were not 
made for our sources, and therefore we cannot apply any background subtraction to our 
high-resolution data. As a result, the continuum of the high-resolution spectra will always 
contain a background component. 

The spectra were returned to their rest frames using the redshifts of the targets obtained 
from the NASA Extragalactic Database (NEDS)@. The continuum flux was calculated from 
the local continuum in the low resolution data. Calculation of the continuum luminosity was 
done with distances determined from z where we assumed a value of 73 km/s/Mpc was used 
for H (Riess et al 2005). 

Line strengths were measured from the high-resolution spectra. Since a significant 
number of our sources only have SH data (9.9/zm < \sh < 19.6/xm), we will focus on studying 
the narrow-line Baldwin Effect in strong emission lines in the SH wavelength range. More 
specifically these lines and their corresponding ionization potential are: [SIV] 10.501/im 
(47.22eV), [Nell] 12.814/xm (40.96eV), [NeV] 14.320/im (126.25eV), and [Nelll] 15.555/im 
(63.42eV). Line strengths were determined by fitting a Gaussian function to the feature at 
the central wavelength. The local continuum was subtracted off with a linear fit. 

Equivalent widths were measured by dividing the line strength measured in high-resolution 
by the continuum strength measured in low-resolution. This was done because the low- 
resolution data do not contain the background flux present in the high resolution data and 
thus they more accurately reflect the actual continuum flux. In Table 1, we present line 
strengths, continuum values, and EW for all our objects as well as redshifts and luminosity 
distances. 



2 This research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated by the 
Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics 
and Space Administration. 



Table 1. Various measured properties for the galaxies in our sample. D L is given in Mpc, fluxes are in 
(10~ 21 W cm~ 2 ), EW are measured in (/im), and continuum values are measured in (10 -21 W cm~ 2 fim^ 1 ) . 
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0.114 


482 
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0.0500 


85.4 


N/A 


N/A 


N/A 
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0.0233 
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3.02e3 
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0.0573 


747 
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0.0291 
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33.7 


33.3 


0.0751 


444 


12.5 


0.0239 


525 


31.5 


0.0576 


547 


35.4 


0.0606 


584 


445 


NGC4151 


0.00332 


14.0 


92.7 


0.0202 


4.58c3 


134 


0.0328 


4.10e3 


77.1 


0.0195 


3.94e3 


204 


0.0510 
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0.0913 


889 
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98.0 
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953 


1.04e3 
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9.40 
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33.5 
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32.1 
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0.0693 


112 
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25.3 
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134 
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133 


14.7 
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150 


8.21 
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22.4 


0.124 
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<1.49 


<0. 00268 
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5.19 
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638 
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625 


4.43 
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486 
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0.00618 


26.1 


37.8 


0.0219 
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0.0455 


887 


25.4 
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855 


55.0 


0.0615 


895 


561 
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3. Results 

In Figure 1, we present the EW of the four lines as a function of the nearby mid-IR 
monochromatic continuum luminosity. A narrow-line Baldwin Effect, where EW declines as 
luminosity increases, appears to manifest itself in all four lines across six orders of magnitude 
in luminosity. As mentioned earlier, our sample contains both type I (quasars and Seyfert 
Is) and type II (mostly Seyfert 2s) AGN and we include both types in our analysis. In 
all figures we indicate type 1 objects with filled red circles and type 2 with open green 
circles. We do not see any evidence that either type 1 or type 2 AGN deviate from the 
trend and while type 2 AGN are generally in the lower mid-IR luminosity range, they do 
not skew the overall Baldwin effect. The blue open triangles denote galaxies with significant 
circumnuclear starburst contribution, classified in a manner described below. We exclude 
them from further correlation analysis. 

In Fig. 1 we also include the error bars to our data. The uncertainties in the continuum 
luminosity were generated by calculating the variance in the local continua. No assumptions 
were made about uncertainties in the redshift measurements and therefore our errors are 
likely lower limits. The flux density uncertainty was a combination of variance in the local 
high resolution continuum as well as error in the profile fitting. Only lines detected with 
confidence greater than 3(Jfi ux , were included in our analysis. The error in EW was the 
propagation of the noise both in low resolution (J con ti nuum , and in high resolution <Jfi ux . For 
nearly all our data, the uncertainty was <10% indicating good quality of our measurements. 
The reason why the error bars appear so small on the plot has much to do with the large 
dynamic range in luminosities and EW measured in our sample. 

We used the conventional definition of the Baldwin Effect expressed as 

W x = aL{ (1) 

where W\ indicates the EW and /3 the slope of the anti-correlation. We performed a 
least squares fit to the data in logarithmic space and to test the linear correlation of the log 
of EW versus log of luminosity, a Spearman rank correlation was employed. The best fits to 
the data are also plotted in Figure 1. 

In three lines, [SIV], [Nell], and [Nelll] we find a greater than 3a significance in the anti- 
correlation. For [SIV], the Spearman rank correlation strength is -0.60 with a null-hypothesis 
value of 4.16 x 10~ 5 . The scatter is comparable to the CIV Baldwin effect (Kinney et al. 
1990) and the least-squares slope is —0.29 ±0.05; steeper than the CIV Baldwin effect which 
is near -0.17 (Kinney et al 1990, Wilkes et al 1999, Osmer and Shields 1999, and Laor et al 
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0.0311 
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16.4 
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0.0163 
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2.07e3 


207 


0.0883 
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15.0 
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2.02o3 
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2.15e3 
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840 
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Fig. 1 — Plots of the EW of [SIV], [Nell], [NeV], and [Nelll] as a function of the continuum 
monochromatic luminosity next to the line. The green open circles denote Seyfert 2 galaxies 
while the red filled circles indicate the Seyfert 1 galaxies. Starburst galaxies are marked in 
blue open triangles. The data are overplotted with a lines representing the least-squares fits. 
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1995). The [Nell] anti-correlation also manifests itself with a Spearman rank value of -0.48, 
significance of 8.88 x 10~ 4 , and a slope of —0.25 ± 0.06. For [Nelll], the correlation is -0.46, 
1.17 x 10~ 3 significance, and a slope of -0.22 ± 0.06. 

The significance of these correlations were tested in two different ways. The first was 
by calculating the number of deviations from the null-hypothesis value for the Spearman 
rank coefficient. The second was by running Monte Carlo simulations where we randomly 
assigned our EW values to our continuum measurements. This was done 100,000 times for 
each line and the resulting correlations were noted. In all cases where we claim to have a 
greater than > 3a significance, the Monte Carlo simulations provided a second, independent 
verification. 

Unlike the first three lines, when we examine the [NeV] line we find that the correlation 
strength is just -0.41 (5.34 x 10~ 3 significance) with a slope of —0.19 ±0.06. The 2.5cr signif- 
icance and low correlation value are considered marginal. As we will discuss in the following 
sections, this was rather unexpected since [NeV], due to its high excitation potential, is tell- 
tale sign of AGN activity. Since the Baldwin Effect is, at least in the optical/UV, directly 
coupled with the AGN activity one would expect that it would be more prominent in strong 
mid-IR AGN lines as well. 



3.1. Eddington Luminosities 

In order to examine the possibility of the Eddington ratio (L/LEdd) as a driver for the 
mid-IR narrow line Baldwin effect, we searched the literature and calculated L/L Edd for 29 
of our sources that had measured bolometric luminosities and black hole masses. 

Our analysis revealed no strong relationship between Eddington Luminosity and EWs 
for mid-IR lines. The Spearman rank correlation values we find are -0.21 for [SIV], -0.45 for 
[Nell], -0.36 for [NeV], and -0.36 for [Nelll]. Their corresponding null-hypothesis values are 
0.320, 0.032, 0.070, and 0.063. The low correlation values coupled with the fact that none 
of relations have > 3<r significance indicates that the Eddington ratio is not a driving factor 
for the weakening of lines in the narrow line region. 

3.2. Slopes of the Baldwin Effect 

Previous studies of the broad line Baldwin effect have revealed that as the ionization 
potential increases, the steepness of the anti-correlation also increases (Osmer & Shields 
1999, Wu et al. 1983; Kinney et al. 1987; Kinney, Rivolo, & Koratkar 1990; Badwin et 
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al. 1989; Zheng et al. 1997; Dietrich et al 2002). This was first noticed in the comparisons 
between low ionization lines such as Mg II and Lya and high ionization lines such as C IV 
and O VI. 

To examine the relationship between slope and ionization potential in our study, we 
have compiled a list of (3 slopes from the literature (Kinney, Rivolo, & Koratkar 1990; Zheng, 
Kriss, & Davidsen 1995; Dietrich et al. 2002, Croom et al. 2002) and also included data 
from our sample. The result is given in Fig. 2 and demonstrates that the values from our 
narrow- line study of AGN in the mid-IR do not appear confirm earlier assessment that the 
broad-line slope of the Baldwin effect becomes steeper when one considers lines of increasing 
ionization potential. Further confirmation of this finding is given by Honig et al. (2008), 
who presented anti-correlations between mid-IR EW's and 2-10 keV luminosity. They found 
that the steepness of their anti-correlations were also independent of ionization potentials. 

Comparing to narrow-line optical data, the results are inconclusive. Croom et al. (2002) 
only reported strong anti-correlations for two lines, [NeV] and [Oil] (plotted as blue, open 
diamonds in Fig. 2), but felt that [Oil] might be contaminated by host galaxy emission. 
Therefore, there are not enough data points to draw definitive conclusions. They suspected 
that optical, narrow line slopes steepen with ionization potential and, if so, then the physical 
processes driving mid-IR narrow lines may be quite different from their optical counterparts. 

4. Discussion 

In the original Baldwin effect, the EW of the lines was anti-correlated with the strength 
of the UV continuum emission, and the latter was used as tracer of the AGN intensity. This 
technique is not as direct in the mid-IR because the local continuum may not reflect AGN 
intensity and could introduce some scatter in the correlations we present in Fig. 1. It is well 
known that in the infrared, the spectral signatures of an AGN can be severely blended by 
emission originating from circumnuclear starbursts (see Laurent et al. 2000, Armus et al. 
2007 and references therein). This is only partially due to limited spatial resolution of infrared 
telescopes and focal plane arrays. More importantly it is the intervening dust which fully 
reprocesses the intrinsic radiation from the sources - both massive stars and/or an accretion 
disk - and re-emits it in the infrared. This may lead to a difference between the optical 
and mid-IR classification of a source. To address this issue a number of diagnostic methods 
have been developed of the years in order to obtain a robust AGN/starburst classification 
of mid-IR spectra (see Lutz et al. 1998, Laurent et al. 2000, Armus et al. 2007, Spoon et 
al. 2007, Charmandaris 2008, Nardini et al. 2008). 
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Fig. 2. — Plot of ionization potential versus anti-correlation slope, (3. The red filled diamonds 
denote narrow-line data from this paper and have their associated error bars. The blue open 
diamonds are narrow line data from Croom et al. (2002), and black squares are points taken 
from various sources and represent broad-line data. 
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It is important for our analysis to explore whether the observed anti-correlation we find 
is indeed related to the strength of the AGN. Ideally one could address this by examining 
the correlation between the line EW and the strength of the X-ray emission (see Hoenig et 
al. 2008). However, these data are not available for our sample so instead, one can examine 
how strongly the AGN contributes to the mid-IR continuum emission near the lines versus 
how much is contaminated by circumnuclear star formation. 

It is widely believed that [NeV] 14.32 /xm, due to its high ionization potential, originates 
in the NLR (see Gorjian et al. 2007 and references there in). So one would have to examine 
whether most of the emission from [Nell], [Nelll], and [SIV] in AGN also come from the 
same region. Gorjian et al (2007) found that both the [NeV] 14.32 /iin and [Nelll] 15.6//m 
are strongly correlated and deduced that the two must be produced in the same region. 
Testing the correlation between [NeV] and [Nelll] in our sample confirms this results at 
> 5o- confidence level. We also examined whether [SIV] and [Nell] correlate with [NeV] 
and found that as with [Nelll], both possessed > 5a relationship with very little scatter. 
Altogether, this suggests that the four lines investigated in this paper are likely to arise from 
the same region. 

To estimate the contribution of massive star formation in the circumnuclear regions of 
the AGN sampled by the IRS slits, we used two diagnostics which have been proposed by 
Genzel et al. (1998) and have also been applied by Armus et al. (2007) in a study of local 
ultralumninous infrared galaxies. The first diagnostic was to use the [NeV] /[Nell] ratio. This 
ratio is useful because it is very difficult for stellar sources to produce a significant number 
of photons capable of ionizing [NeV]. As a result a high [NeV] /[Nell] line ratio indicates the 
presence an AGN which is dominant in the mid-IR (Lutz et al. 1998). 

A second diagnostic is to use the EW of the 6.2/im emission feature. When this feature 
is strong, it indicates the presence of Polycyclic Aromatic Hysdrocarbons (PAHs) which are 
mainly produced in photo-dissociation regions (PDRs) when they are excited by the adjacent 
star forming regions. Therefore, in Figure 3 we plot the EW(6.2/im) versus [NeV] /[Nell] 
similar to Figure 5 in Armus et al. (2007). From Figure 3 we see that there are eight 
galaxies (NGC 7469, IC 5135, NGC 5135, Mrkl066, NGC 1275, NGC 2110, NGC 2273, and 
NGC 7213) that have both a low [NeV] /[Nell] value and strong PAH emission. As mentioned 
before, we flagged them as starburst galaxies and did not use them in our correlation analysis. 

Even after removing the most flagrant starburst contaminants, it is evident that Fig. 1 
still possesses a scatter. This is largely due to the fact that star formation activity is not 
a binary phenomenon where it is either overwhelmingly dominant or completely dormant. 
The Armus et al. (2007) and Genzel et al. (1998) figures show that the diagnostics used 
are continuous between the two extremes and demonstrate that while there are cases where 
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Fig. 3. — The EW of the 6.2/im PAH feature versus [NeV]/[NeII]. The red diamond indicates 
3C273, a well known quasar, where an AGN dominates its mid-IR spectrum (Hao et al. 2005). 
The blue open triangles denote objects with both high PAH emission and a low neon ratio. 
These are the strongest candidate starburst galaxies and were flagged as such in our sample. 
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star formation or the AGN activity dominate, in many galaxies the signature of both in the 
mid-IR is of similar strength. Therefore, even though we removed the galaxies where we 
were certain starburst activity was dominating, the remaining classified as AGN still possess 
some scatter in our anti-correlations. This is very likely one of the reasons why our [NeV] 
result is not strong. 

It is evident from the figures that the two diagnostics used are continuous between 
the two extremes and demonstrate that while there are cases where star formation or the 
AGN activity dominate, in many galaxies the signature of both in the mid-IR is of similar 
strength. Therefore, even though we removed the galaxies where we were certain starburst 
activity was dominating, the remaining classified as AGN still possess some scatter in our 
anti-correlations. This is very likely one of the reasons why our [NeV] result is not strong. 

To reduce the possible contribution of star formation to the observed correlations further 
yet, we removed the local continuum from the analysis altogether. Laurent et al. (2000) 
and Nardini et al. (2008) point out that the 5.5/im continuum is dominated by emission 
from dust located near the AGN torus which heated to near sublimation temperatures, with 
minimal contribution from star formation activity or stellar photospheric emission. Since 
the flux of a mid-IR line is independent of local continuum levels, the ratio of the line flux 
to the 5.5/im continuum versus 5.5/zm luminosity should mitigate continuum contamination 
altogether. 

The result of this experiment is clear. The overall ant i- correlations are better than those 
of Fig. 1. All four lines display > 3a significance and [Nell] and [Nelll] actually have > 4<r 
significance. The correlation values are -0.58 for [SIV], -0.65 for [Nell], -0.57 for [NeV], and 
-0.62 for [Nelll], with null hypothesis values of 1.03 x 10~ 4 , 1.29 x 10~ 6 , 6.41 x 10" 5 , and 
3.16 x 10~ 6 respectively. Using this method, the slope values are —0.32 ± 0.06 for [SIV], 
-0.36±0.05 for [Nell], -0.31 ±0.06 for [NeV], and -0.37±0.06 for [Nelll]. As expected the 
starburst galaxies show some of the strongest [Nell] emission in our sample and are obvious 
outliers. The figure further suggests that the anti-correlation is present in all four lines and 
is driven by the central AGN. The relationship probably only weakened by star formation 
in the host galaxy. These results are more meaningful than the ones from Fig. 1 principally 
because of the inclusion of the [NeV] line which, as mentioned before, we expect most to 
exhibit a Baldwin effect due to its close relation to the central AGN. Therefore, its inclusion 
in this analysis points to the AGN as the central driver of this effect. The analysis used in 
Fig. 4 is shown to be a robust way to measure the decrease in line strength with increasing 
AGN power and may be a better diagnostic tool than the traditional analysis used in Fig. 1 
for mid-IR lines. 

Altogether, our data support the possibility that the driver of the narrow-line Baldwin 
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Fig. 4. — Plot of line flux divided by 5.5/xm flux versus 5.5/zm Luminosity for [S1V], [Nell], 
[NeV], and [Nelll]. The green open circles denote Seyfert 2 galaxies while the red filled circles 
denote Seyfert 1 galaxies. Starburst galaxies are plotted as open blue triangles. Overplotted 
are lines representing the best least-squares fits. 
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effect is a dynamic covering factor of the narrow line region. If this factor changes as a 
function of AGN luminosity where more luminous AGN tend to drive the NLR outward, 
then it would explain the decrease in relative lines strength with increasing AGN power seen 
in Fig. 1 and 4. 
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